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We reappraise and clarify the contradictory statements found in the literature concerning the time- 
of-arrival operator introduced by Aharonov and Bohm in Phys. Rev. 122, 1649 (1961). We use 
Naimark's dilation theorem to reproduce the generalized decomposition of unity (or POVM) from 
any self-adjoint extension of the operator, emphasizing a natural one, which arises from the analogy 
with the momentum operator on the half-line. General time operators are set within a unifying 
perspective. It is shown that they are not in general related to the time of arrival, even though they 
may have the same form. 
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I. INTRODUCTION 

Even though there is no question that experimental- 
ists measure distributions of time of arrival for quantum 
systems, there has been a long debate about the capa- 
bilities of standard quantum mechanics to address the 
very concept of time as an observable. In a nutshell, 
the problem arises because we expect observables to be 
represented in the quantum formalism by self-adjoint op- 
erators, whereas an old theorem of Pauli |Q] states that 
for a semibounded self-adjoint operator H no conjugate 
self-adjoint operator T can exist, i.e., no operator that is 
self-adjoint and satisfies the canonical commutation re- 
lation [H, T] = ih over a dense domain. In other words, 
there is no self-adjoint time operator if the Hamiltonian 
is bounded from below, as we normally expect it to be. 
However the theorem has not discouraged theorists from 
attempting to fit such an immediate classical concept into 
the standard framework (see for a general discussion) . 
In relatively recent times a rather satisfactory answer has 
been given for the arrival time of freely moving states 
in one dimension, from the point of view of generalized 
spectral functions, or positive-operator- valued measures 
(POVM) H [| . In short, the relevant statement is that 
sclf-adjointncss is not strictly necessary for the formula- 
tion of probabilities and the reality of measurements of 
observables, and that, in fact, any given observable is 
uniquely characterized (in a specific sense) by the prob- 
ability distributions of measurement results in the dif- 
ferent states accessible to the system (for a much more 
developed presentation, refer to QJ|]). 

This result has sparked a renewed interest in the anal- 
ysis of old and new proposals for time operators in gen- 
eral, and in particular for the time-of-arrival operator 
Tab (see below) introduced by Aharonov and Bohm 0, 
as in §§|§. ' 

However, a number of basic issues are still in dispute, 
or in need of clarification. In particular, contradictory 
statements may be found about the domain of the time 



operator Tab or a bout its eigenfunctions ]l6||i"4| ]; 

and the deficiency indices, which arc important to deter- 
mine its self-adjoint extensions, appear to be, following 
different authors, either (0,1) or (0,2) 00, but only 
one of the possibilities may be correct. Even though sev- 
eral self-adjoint "variants" of Tab have been proposed, 
|ll|]l6| , |l5| none of them is a self-adjoint extension of Tab 
in a proper sense. The framework to consider these ex- 
tensions is Naimark's dilation theorem PJ^,p^|, but an 
explicit simple construction of a self-adjoint extension is 
lacking (For a rather abstract presentation of the the- 
orem applied to the time-of-arrival problem, under the 
guise of "screen observables", see also [f|]|,||])- In 
addition, the form of the time operators in different rep- 
resentations may also be, and has been, misleading: the 
same operator has been presented in different forms, but 
different operators have also been written in the same 
form. 

It appears that in order to make further progress in 
this field without being trapped by any of the numerous 
pitfalls mentioned, a detailed analysis of these topics is 
required, and their clarification is the main objective of 
this paper. Even though some of the arguments have 
to be rather technical, our aim is to provide an account 
readable by a general physicist. We shall in particular in- 
troduce a description of the extension theorems as they 
apply to the problem at hand, and develop a general 
analogy, sketched by Holevo Q , between the momentum 
operator in the half line and the time operator. We hope 
that it will provide a much better understanding of the 
relevant issues. It will also allow to join the strands given 
by i) the initial proposal of the time-of-arrival operator 
of Aharonov and Bohm 0, its modifications or variants 
|Ti]|l|-[17| ; ii) the POVM idea as applied to quantum 
physics and to the time operator in particular iii) 
the axiomatic approach of Kijowski to obtain good prob- 
ability distributions for the time of arrival Jll]|l^] ; and iv) 
some apparently formal constructions in AUcock's semi- 
nal paper S. 
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As is well known, the momentum operator in the hair- 
line is is not self-adjoint. Moreover, it is maximally sym- 
metric, because its deficiency indices are (1,0), and it 
admits no self-adjoint extension on the Hilbert space of 
functions defined on the half-line (for a fast and good 
review of deficiency indices and the Neumann theory of 
extensions, with particular emphasis on the physical ex- 
amples of momentum operator and Laplacian, see J2j| - 
we will anyhow give a self-contained description in sec- 
tion II). On the other hand, there is a natural exten- 
sion, namely, the momentum operator on the full line. 
It should be noticed that this is most definitely not the 
unique extension, but just one of the infinitely many pos- 
sible. Nonetheless, its "naturalness" is well justified from 
a physical point of view. The question might then be 
posed as to how are we to recover any information about 
the distribution of probability of momentum in the half- 
line from the full line extension. This is actually the 
content of Naimark's dilation theorem, together with the 
theorem that asserts the uniqueness of the POVM asso- 
ciated with a maximally symmetric operator: that the 
unique distribution of probability associated with the 
momentum operator on the half-line is recovered from 
the spectral decomposition of the self-adjoint operator 
momentum on the full line, by simple projection. Equiv- 
alently stated, we can reproduce the physical content of 
the momentum operator on the half-line from the stan- 
dard, well-known techniques that apply to self-adjoint 
operators. 

We will show that the different alternatives that have 
been proposed to make sense of the time operator of 
Aharonov and Bohm are actually consequences of this 
general result, that we have phrased in terms of momen- 
tum operators, but which can be extended to other cases 
as well. The generality of the result also explains the 
coincidence of different constructions, in particular the 
axiomatic approach of Kijowski [|ll|Jl2| and the attempts 
to create ex novo a self-adjoint operator ( jl^] is a rele- 
vant example). 

We organize the paper as follows: in the next section 
we give a complete and very basic description of the mo- 
mentum operator on the half-line and the problems that 
appear there. We also give the natural extension and 
describe more fully Naimark's dilation theorem in this 
context. In section III we establish the parallel with the 
time-of- arrival operator of Aharonov and Bohm, to be 
followed in the next section with the general case of a 
time operator. Section V is devoted to discussing and 
restating the different results known in the literature in 
the terms used in the previous sections. 

II. THE MOMENTUM OPERATOR 

Consider the momentum operator p = —ihd x defined 
on a dense domain of the Hilbert space of square in- 
tegrable functions on the half-line, 7i> = L 2 (R + ,d:r). 



More precisely, it is defined on D(p), the subspacc 
of square integrable absolutely continuous functions ip 
whose derivative is also square integrable, and such that 
ip(0) = 0. Clearly, this operator is symmetric on its do- 
main. Let us write down a set of generalized (weak) 
eigenfunctions that provide us with a resolution of the 
identity. The adequate set, which does not belong to the 
Hilbert space and is parametrized by real p, is given by 
ip p (x) — eKp(ipx/h)/V2nh, as is only to be expected. 
They form a complete basis, since dpip p (x')ip p (x) — 
5{x' — x). However, they are not orthogonal (in the gen- 
eralized sense): 

r°° i i i 

/ dxip p >(x)ip p (x) = -6{p-p) + — P 7 , 

Jo 2 2vr p-p' 

where P stands for principal part. This is evidently dif- 
ferent from the usual constructions for self-adjoint oper- 
ators. 

Let us now pass on to the position operator x, defined 
on those square integrable functions on the half-line, ip g 
H>, such that J °° dx x 2 \tjj(x)\ 2 < +oo. This operator is 
self-adjoint on its domain. Additionally, it is true that 
[x,p] = ih on a dense domain. However, in the case 
at hand the position operator is bounded from below, 
and Pauli's theorem therefore applies. It follows that p 
is not a self-adjoint operator and admits no self-adjoint 
extension in H>. 

We shall now rephrase the statements of the previous 
two paragraphs in a different way. The domain of the op- 
erator adjoint to p, i.e., p\ according to von Neumann's 
formula, is D(p) ® N(i) N(-i), where N(±i) are the 
spaces of eigenvectors of p* with eigenvalues ±i respec- 
tively. This comes about because, even though p, be- 
ing symmetric, has no imaginary eigenvalues, its adjoint 
does have eigenvectors with imaginary eigenvalues. As a 
matter of fact, the functions exp(i\x/h), with A a com- 
plex number with positive imaginary part, are all of them 
eigenvectors of jft of eigenvalue A. Notice that the dimen- 
sion of N(X), the space of eigenvectors with eigenvalue A, 
is the same for all A with positive imaginary part. This 
is called the first deficiency index of the operator p, while 
the second one is the dimension of the space of eigenvec- 
tors for any given eigenvalue of negative imaginary part, 
in our case, 0. That is, the deficiency indices of p are 
(1,0). But self-adjoint operators have deficiency indices 
(0,0), as do essentially sclf-adjoint operators, which, al- 
though not self-adjoint, have a unique self-adjoint exten- 
sion, namely, their closure. If we were to build a self- 
adjoint extension of p we would need to include somehow 
those elements of D{p*) that are not in D(p), but in such 
a way that the imaginary parts compensate. When the 
deficiency indices are equal, this is achieved by the use of 
a unitary transformation from N(X) to N(X) (this is the 
von Neumann theory of self-adjoint extensions, whereby 
each self-adjoint extension in the Hilbert space of defi- 
nition of an operator over a dense domain, with equal 
defect indices, is given by a unitary transformation be- 
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tween the deficiency subspaces) . As the deficiency indices 
in our case are different, it is not possible to extend p to 
include an action over the whole of the domain of D(p s ) 
such that the imaginary parts compensate. Since p is 
a closed operator, symmetric over its domain, which is 
dense in 7i>, and its deficiency indices are unequal, one 
of them being 0, we say that it is a maximally symmet- 
ric operator, and the previous result is that it admits no 
self-adjoint extension over 7i>. 

Retaking the complete set of generalized eigenfunc- 
tions ip p (x), we can construct a positive-operator- valued 
measure (POVM) F. This is a map from intervals in the 
real line (a c-algebra of subsets of a nonempty set) to 
the positive operators over a Hilbert space that satisfies 
three axioms that we will illustrate with the example at 
hand. The map for our case is given by 



■\F([a,bM 



dp 



dx 



dy (j>(x)ipp(x)ipp{y)ip{y) . 



over all <f>, ip S H > . This map sends intervals of the 
real line to positive operators acting on 7i > (first axiom) , 
which add together when the intervals are disjoint (sec- 
ond axiom), and which add to the identity operator when 
summed over the real line, because of the completeness 
proved above (third axiom). This differs from the usual 
decompositions for self-adjoint operators in that it does 
not fulfill the property that the positive operators be pro- 
jectors, i.e., in our case F([a,b]) 2 ^ F([a,b}), because of 
the lack of orthogonality shown above. 

Now we have the POVM associated with p, we can 
reconstruct the operator and the probability distribution 
of its values. The action of the operator is given by 



(?¥>)(%) = / d P dypip p (x)ip p (y)<p(y) . 

J -oo JO 

The probability distribution for its possible values over a 
state ip is 



n*(p) = 



dx ip p (x)ip(x) 



Notice that this indeed satisfies all the requirements for 
it to be a probability. 

Let us see in more detail the relationship between the 
probability density and the POVM. We have written 
above the expression for the action of F([a, b]). Sym- 
bolically, we can write F(dp), and its expectation value 
on a state ip is 



(ip\F(dp)ip) = dx dyip{x) 
Jo Jo 
= U i ,(p)dp. 



Jp(x-y)/h 



2-kTi 



-V>(y) dp 



The expression just written defines in fact (4>\F \dp)ip) , 
and therefore the whole POVM. Furthermore, the do- 
main of the operator associated with the POVM is also 



defined by IL^ (p): by integration by parts it can readily 
be seen that the domain of the operator p is precisely the 
set of states ip for which the second moment of the prob- 
ability distribution H^(p) is finite, and there it coincides 
with ||jh/;|| 2 . Therefore, the probability density written 
above defines the POVM and the operator with it asso- 
ciated. Notice also that the probability density need not 
be a continuous function. The only requirement is that 
ILy,(p)dp be a good measure on the real line. 

Another important property of the POVM is covari- 
ance under displacements of momenta, which reflects the 
commutation relation [x,p] — ifi. Namely, it is readily 
computed that, for all real q, 



lq *f h F{[a,b])e- lq± / n ij) 



\F([a + q,b + q})iP) 



In terms of the probability density, the statement is that 



where ip q 



iqx jh 



ip is the shifted state. 



The probability density, written above for the case of 
pure states, can be easily generalized to mixed states. 
Let p(y, x) be the matrix elements of the density matrix 
p in the position representation. The probability density 
associated with this density matrix is then 



d.r 



e ip(x-y)/h 

dy r-r — p(y,x) 

o 2ttTi 



In what follows we will use only pure states in the dis- 
cussion, bearing in mind that the generalization to mixed 
states is straightforward. 

It would seem that something akin to the spectral the- 
orem of self-adjoint operators has been achieved, and this 
is indeed the case. The difference, however, lies in that 
the expectation value of higher order powers of p on a 
state ip, (ip,p n ip), does not necessarily coincide with the 
corresponding moment of the distribution ILy, (p), i.e., 
J^ y2 dpp n H^(p), for n > 3. For the case at hand, for 
instance, this happens already for n — 3 if ^'(0) 7^ 0. As 
a matter of fact, in such a situation the expectation value 
for p 3 has an imaginary part. Nonetheless, the measures 
of momenta can be readily associated with ILy,(p), which 
carries the relevant physical information. 

Even so, we would normally like to understand the con- 
structions above in terms of the more usual recipes for 
self-adjoint operators. We shall now make essential use 
of the uniqueness theorem for POVMs of maximally sym- 
metric operators: given a maximally symmetric operator 
A over a Hilbert space H, there is a unique POVM, Fa 
(unique up to isomorphisms), such that its first operator 
moment coincides with the operator, and that the set of 
states over which the second moment exists is precisely 
the domain of A H^2|. This means that if we are to 
construct by whichever means a POVM such that it ful- 
fills these conditions, we will be obtaining again the same 
POVM. 
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Moreover, Naimark's dilation theorem tells us that any 
POVM associated with a symmetric operator defined on 
a dense subset of H can be constructed from a self-adjoint 
extension of the operator to a larger space 7i, as follows: 
let E be the projection valued measure of the self-adjoint 
extension (i.e., a POVM that satisfies the further require- 
ment that E([a, b]) 2 — E([a, &])), and P the projection 
operator from Ti to TL. Then F([a,b]) := PE([a,b]) is a 
POVM associated with the symmetric operator. 

As a consequence, if we are to build a self-adjoint ex- 
tension of p in a larger space, we reproduce the unique 
POVM and the whole of the physical content of the oper- 
ator from the usual analysis for the self- adjoint extension 
and a projection. Notice however that the possible ex- 
tensions are infinite. Not so the POVM, and that makes 
the freedom of choice of extension even more interesting. 

Back to the case of p — —ihd x defined on the half-line, 
we see that there is a simple possibility: to extend the 
operator to the full line, i.e., —ihd x defined on a dense 
subset of H = L 2 (R, dx). Naturally enough, the action of 
this operator on the elements of D{p) is the same as that 
of p, so it is an extension. This is not the only self-adjoint 
extension on the whole real line, of course. As a matter of 
fact, the deficiency indices for the direct sum of the mo- 
mentum operators on the positive and negative half-lines 
are (1,1), thus signalling that a one dimensional contin- 
uum of alternative self-adjoint extensions exist. They dif- 
fer by the presence of a jump function located on x — 0. 
However, most natural is the momentum on the full line, 
defined over absolutely continuous functions, i.e. with no 
jump on x — 0. This is a self-adjoint operator, for which 
the standard spectral analysis is applicable. 

More concretely, the projection valued measure E is 
given by the expression 

rb poo roc 

{4>\E{[a,b])f) = / dp dx dy^)%{x)W)^{v) 

J a J — oo J — oo 

with ip p (x) = exp(ipx/K)/ VZnh, as before, but now de- 
fined over the whole real line. The projection P in our 
case is simply (Pip)(x) = <d(x)ip(x), with O being Heav- 
iside's step function. 

The probability distribution for the momentum oper- 
ator over the whole real line for a state ip is of course the 
modulus squared of the wave function in the momentum 
representation, and its restriction to states that belong 
to D{p) is none other than the probability distribution 
associated with the POVM. 

We see then that the reason we could not do the stan- 
dard analysis for the momentum operator on the half-line 
is that we are being, in a way, far too restrictive in the 
behaviour near x — of the states on which it can act. 
There is a reminder of the full line, seen for instance in 
the principal part that forbids orthogonality, or in the 
fact that the higher moments of the probability distribu- 
tion do not in general agree with the expectation value 
of powers of the restricted momentum operator. If our 



states were such that the function ip and all of its deriva- 
tives were zero at x = 0, there would be no problem 
with the higher moments, and we would get no imagi- 
nary part for the powers of p over such states. However, 
such a strong restriction on the allowable wavefunctions 
would cut out many physically sensible cases. 

Let us now write slightly more formally the extension 
procedure we have performed to reobtain and interpret 
the probability density: we have started with a maxi- 
mally symmetric operator p = —ifid x acting on the dense 
domain D(p) C 7i> = L 2 (R + ,dx), with deficiency in- 
dices (1,0). We have then considered an extension in 
L 2 (R, dec) , making use of the isomorphism 

L 2 (R, dar) = L 2 (TL+,dx) ® L 2 (Br, da) . 

The extension has been the natural one, i.e., —ihd x on 
the full line (although, as we repeatedly stated, this is not 
the only possible extension; another easy choice would 
have been some extension of p® (—p) , for instance). The 
standard spectral analysis for this operator produces the 
probability density, which, when restricted to 7i>, gives 
n^(p), the probability density out of which the POVM 
for the operator we started from can be built. 

III. THE TIME-OF- ARRIVAL OPERATOR 

Within the long-running discussion on the concept of 
time in quantum mechanics, and in particular with re- 
spect to the status of the time-energy uncertainty rela- 
tions, Aharonov and Bohm introduced in an important 
paper Q, among other questions, "a clock" to measure 
time from the position and momentum of a freely moving 
test particle. The corresponding operator was obtained 
by a simple symmctrization of the classical expression 
' my/py, where y &ndp y are, respectively, the position and 
momentum of the test particle. By the same token, the 
operator obtained by symmetrizing the classical expres- 
sion for the arrival time at x = of a freely moving par- 
ticle having position x and momentum p, t = —mx/p =, 
is given by |@ 

Tab ■= -y [xp +p x) , 

Note the minus sign in comparison to the clock time in 
|7| . In spite of the somewhat subtle difference (in concept 
and sign) with the original time operator introduced in 
we shall refer to this operator as the Aharonov-Bohm 
(time-of-arrival) operator. Without regard for topologi- 
cal considerations, it is clear that it has the correct com- 
mutation relation with the free particle Hamiltonian on 
the line, Ho = p 2 /2m. It is thus a good candidate for a 
time operator, with the plausible physical interpretation, 
given by the correspondence rule, that it is related to 
the time of arrival (notice that other quantization rules 
might possibly give a different result, although this is 
the operator obtained not just by the symmetrization 
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rule, but also by applying Weyl, Rivier, or Born- Jordan 
quantizations f24|| ). However, it also follows, from Pauli's 
theorem, that it cannot be a self-adjoint operator acting 
on a dense space. 

Let us then examine first the question of the domain 
of Tab- In order to do that, it is useful at this point to 
consider the Hilbert space of the free particle in the mo- 
mentum representation, 7i p := L 2 (R, dp) (on this space 
we know how p _1 acts, thanks to the spectral theorem). 
Formally, we then obtain 



Tab 



ihm ( 1 



P 



2d_ 

p dp 



There are other alternative expressions, such as 
—ihmp~ 1 / 2 dpp~ 1 / 2 , which would be valid for p > 0, or 
for p < by analytic continuation, see e.g. [Q. At any 
rate, Tab understood as a differential operator presents 
a singular point at p = 0, and this is the source of all the 
difficulties that have appeared in the literature. 

The differential operator written above can only be 
applied to absolutely continuous functions, but there are 
further requirements. One such is that TabV' belongs 
to Hp, i.e., that it be square integrable. This poses 
a restriction due to the singularity of the operator at 
p = 0. On computation, one finds that the singular- 
ity is avoided if the function ip has one of the following 
possible behaviours close to p — > 0: either ip(p) ~ p 1 / 2 , 
or ip(p)/p 3 ^ 2 — > 0. However, this is not enough to fix 
the domain of the operator, as Paul noticed long ago 
jjj. Given that, at least formally, Tab is symmetric, 
this should also be a requirement on its domain. In- 
tegration by parts, and demanding that {ip\Tab^p) be 
equal to (Tab^pIp) for all ip in the domain of Tab, leads 
us to exclude the first possibility, thus defining the do- 
main of Tab, D{Tab), as the set of absolutely continuous 
square integrable functions of p on the real line, such that 
ip{p)/p 3/2 -> as p -> and \\f A B^\\ 2 is finite. As for 
the alternative expression —ihmp~ 1 / 2 d p p~ 1 / 2 , let us ask 
that (©(p) — p))ip(p) / yj\p\ be absolutely continuous 
for ip to be in its domain. This, together with the fur- 
ther requirement of symmetry, leads us to a domain that 
coincides with that of Tab ■ Since the respective actions 
also coincide over this domain, we see that they are but 
equivalent differential expressions for the operator, once 
the adequate domain is taken into account. 

Tab is closed over this domain, as can be checked by 
,t 



computing [T AB j (there is no doubt that D(Tab) is 
dense in 7i p , so the closure of Tab must coincide with 
■Tab") )• Before that, though, let us examine Tab anc ^ 



its domain. 

In order to apply von Neumann's formula, we have to 
check whether there are states ip in H. p such that for all 
ip G D(Tab) the following expression holds: 



(Tab + P) 



since then ip is an eigenvector of T\ B with eigenvalue i. 
Analogously, we also have to study the case with eigen- 
value — i. By integration by parts, and application of the 
condition that all functions in D(Tab) satisfy, namely, 
that ip(p)/p 3 / 2 tends to zero as p tends to zero, it is 
found that there are two independent eigenvectors with 
eigenvalue i, and none with eigenvalue —i. The relevant 
eigenvectors are 



ip±(p) = 0(±p) V±pe 



-p 2 /2mh 



Notice that, contrary to the expectations of some authors 
|Io|]li| , |l7|| , both of them have to be taken into account: 
there are no requirements of derivabihty or continuity for 
the functions in the domain of the adjoint. 

The deficiency indices are therefore (2, 0), and we have 
a maximally symmetric operator. As we have seen in 
the previous section, this implies that no self-adjoint ex- 
tension can exist, and we should redo for this case the 
analysis performed for the momentum operator on the 
half-line. However, the most convenient way of doing 
that is by passing to the energy representation, as was 
already pointed out by Allcock in his seminal work ||, 
and emphasized again by Kijowski This change of 

representation is useful, from the mathematical point of 
view, because it implements a theorem p2| which states 
that a simple symmetric operator with deficiency indices 
(n, 0) can be decomposed as a direct sum of n operators 
with deficiency indices (1,0). Since each of these is in fact 
isomorphic to the momentum operator on the half-line, 
we will be able to use directly the previous results. 

The change of representation corresponds to the de- 
composition of the Hilbert space H. p into the subspaces 
of positive and negative momentum, i.e., 

L 2 (R, dp) = L 2 (R+ dE) 8 L 2 (R+, dE) 

= H+®H-, 

where the first subspace is that of positive momenta, 
whereas the second corresponds to negative p. The ex- 
plicit isomorphism is given by 



tp±(E) = (m/2E) 1/A ifj(±V2mE) , 



ip(p) 



1/2 



0(p)^ 



P 

2m 



0(-p)V- 



P 

2m 



0. 



where ip € TL P , ip± S TL±, and the isomorphism relates 
ip <-> (ip+,ip-). The factor E~ 1 /^ is due to the change 
in the measure from dp to dE and reciprocally for the 

1 /2 

factor (|p I /to) . The interesting point is that given this 
isomorphism, the time operator of Aharonov and Bohm 
takes the form —ihdE, as is well known. The domain of 
the operator, D{Tab) is sent by the isomorphism into the 
direct sums of square integrable absolutely continuous 
functions in each subspace, such that for each subspace 
we have the restriction that ip a (E)E~ 1 / 2 — > as E — » 
0. But this is exactly the case considered above, since 
the requirement ip a (0) = and the square integrability 
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of ip' a {E) imply the restriction stated before. In other 
words, we have the isomorphism 

Tab = {-ihd E ) © {-iW B ) = f+ © f - , 

where T± are isomorphic to the momentum operator on 
the half-line. 

Therefore, the constructions carried out in the previ- 
ous section can be immediately translated to this situa- 
tion, but taking into account that the energy spectrum 
is degenerate whereas the position spectrum is not. For 
instance, the complete non-orthogonal set of generalized 
eigenf unctions ip p is doubled here into a set with a contin- 
uous parameter t (the notation is intended to be sugges- 
tive) and a discrete one, with values + or — , as follows: 



and 



1 



JEt/H^ 



jEt/h 



In what follows we will not make the explicit distinc- 
tion between the clement of the full Hilbert space and its 
component, if only one is zero. 

These functions transform under the isomorphism to 
give the following expressions: 



$?(p) = e(ap)( 



ap 



1/2 



t ip t/2mh 



2TrmhJ 



It is straightforward to prove completeness, i.e., 
E Q Ho dt (P')^ ] (P) = S (P ~ P')- Alternatively 

E Q JZo dt ( E ')^ ] ( E ) = S ( E - E ')^ where 1 is the 

two by two identity matrix (S(E — E')TL is the identity 
operator on the full Hilbert space Ti + Ti ). 
Nonorthogonality is also a direct translation: 

r dE ^? (E)if>® (E) = ±* Qa , (s{t - t') + ^j^j 



or 



It now behooves us to compute the POVM or, alter- 
natively, the probability distribution for measured values 
of the Tab operator from which it can be readily recov- 
ered. By direct translation, we can write the probability 
density as 



n 



(-0+ ></>-) 



e —iEt/h 

dE^=^+(E) 



dE- 



,-iEt/H 



-4>-{E) 



in the energy representation, or as 

P \V2 



dp ( 



ip t/2mh 



ip(p) 



( - \ 

dp ( I e" 

V ZTTmn I 



ip 2 t/2mh 



in the momentum representation. This is, of course, the 
same as Kijowski's probability density. The essential 
property of covariance under transformations generated 
by the Hamiltonian is also evident and a direct transla- 
tion of the covariance property signalled in the previous 
section. Physically, it means that the probability of ar- 
riving at t for a given state is equal to the probability 
of arriving at t ~ r for the same state evolved a time r. 
This is the reflection on the probability density of the 
canonical commutation relation [H, Tab] = ih. 

The domain of Tab, now defined through the POVM, 
can be characterized as the set of elements of the Hilbert 
space for which J_ oa dt I\^(t)t 2 is finite, and this quan- 
tity defines {Tab^Tab^P) (thus realizing the minimum 
variance demanded by Kijowski and Werner p] , p^| ). 

As before, we would like to understand all these con- 
structions in terms of a generalized self-adjoint exten- 
sion, by using the uniqueness of the POVM associated 
with a maximally symmetric operator, and Naimark's 
theorem. From the structure of the time operator of 
Aharonov and Bohm, namely Tab — {—ihdE)® {—iTiBe) 
on I? (R+,cLE) © L 2 (R+,dE), it follows that a natu- 
ral and simple extension (natural in the sense of follow- 
ing the natural extension in the analogy) is the operator 
{-ihd B ) © {-iW E ) acting on L? (R,dE) © L 2 (R, dE), 
which is obviously self-adjoint. In other words, we have 
introduced negative energies, respecting the twofold de- 
generacy of the initial spectrum. On this space the (dou- 
bled) Fourier transform acts as a unitary transforma- 
tion that provides us with the time representation of the 
states, and the probability density for the time of ar- 
rival over a pure state in this extended space is nothing 
but the modulus squared of the wavefunction in the time 
representation. Notice that the time representation cor- 
responds to the doubled space L 2 (R, dt)@L 2 (R, dt). We 
reobtain the by now usual probability density for time of 
arrival, by restriction to the initial Hilbert space. The re- 
strictions to the initial Hilbert space of the applications 
of the spectral theorem are therefore related to H^p (t) , as 
stated before. 

The problems that the vicinity of p = (alternatively 
E = 0) pose for the analysis of the operator of Aharanov 
and Bohm are thus seen to be due to the (physically im- 
posed) restriction to the space of positive energies. There 
is a reminder of the negative energies in aspects such as 
the nonorthogonality of the set of generalized cigenfunc- 
tions, similarly to what happens in the paradigmatic ex- 
ample of the previous section. 



G 



IV. GENERAL TIME OPERATORS 

In previous sections we have discussed the time of ar- 
rival for free motion but in fact its mathematical form, 
and its POVM, can easily be generalized for arbitrary po- 
tential functions. The important point however is that 
in general the resulting operators and POVMs cannot be 
physically associated with an arrival time. 

Consider a self-adjoint Hamiltonian, bounded from be- 
low, defined on a Hilbert space Tt, such that the spectral 
decomposition leads to the isomorphism 

H = ®ti (®]UL 2 (\ai,oo),diH(E))) , 

where all the measures fii are absolutely continuous with 
respect to the Lebesgue measure. That is to say, the 
spectrum of the Hamiltonian is absolutely continuous, 
and for any given E e R the degeneracy is finite, the 
maximum degeneracy being The Hamiltonian 

is realized in each subspace as the multiplication by the 
variable E. Therefore, a natural (but most definitely not 
unique) candidate for the time operator is again —ihdE 
on each of these subspaces, restricted to functions that 
vanish at eij. Again each component of this time operator 
can be extended to the real line, where it will be self- 
adjoint and admit a spectral decomposition, which, by 
projection, will give us again the POVM associated with 
our candidate operator. The POVM, and its associated 
probability density, will be unique for the given operator 
if this is indeed maximally symmetric. 

Notice moreover that if we do extend one of the Hilbert 
subspaces to be square /x-integrable functions on the real 
line, and demand that the covariance property of the re- 
striction of the time operator to that Hilbert subspace 
be mantained on the real line, then, under fairly general 
conditions, the extension of time operator and Hamilto- 
nian on L 2 (R, dfii(E)) will be unitarily equivalent to the 
canonically conjugate pair —ihdE and multiplication by 
E. Therefore, modulo intervening unitary transforma- 
tions, we have a way of constructing POVMs for time 
operators by Naimark dilations. 

The usefulness or otherwise of the particular time oper- 
ator under consideration will then have more to do with 
its actual properties and measured probability density 
rather than with mathematical difficulties that had rid- 
dled the work of many previous workers in the field. 

To illustrate this point further, let us examine the 
standard example of a particle in a constant field. The 
Hamiltonian is H g = p 2 /2m + mgq, clearly seen to be 
unbounded. Pauli's argument is therefore not applica- 
ble, and a self-adjoint operator canonically conjugate to 
H g does indeed exist, namely, T g = p/mg. In the energy 
representation this operator is of course given by —ihdE- 
It can be ascribed to time of arrival to zero momentum, 
but in no possible way to time of arrival to a specified 
position. This example tells us that an operator having 
the form —ihdE in the energy representation, be it as a 
self-adjoint operator, or associated with a POVM, does 



not necessarily mean that we have obtained a time-of- 
arrival operator, even though a time operator is indeed 
being considered. 

At this point it is worthwhile mentioning the approach 
of Leon et al. [E5|, intended to generalize Tab to the 
interacting case of a potential barrier (under the hypoth- 
esis of asymptotic completeness and absence of bound 
states). They do obtain good time operators, involving 
—ihdE in the basis of outgoing or incoming scattering 
states, and their distributions, understood as POVMs. 
However, since the expectation value of these operators 
is, over a given state, the same as the expected time of ar- 
rival that the corresponding asymptotic states (outgoing 
or incoming respectively) would produce in the case of no 
interaction, it is hard to understand them as bona fide 
time-of- arrival operators at an arbitrary spatial point. 

Finally, it is also interesting to examine the relation of 
the operator Tab with older literature on the "time op- 
erator" by Olkhovski, Recami, and others (^6|. The basic 
idea behind this set of works is to extract a time operator 
from the relation that defines an average "presence time" 
at x = by 



(t) 



J+™dt\j;(x = 0,t)\H 
£~<ft|#s = 0,t)| 2 



We shall discuss as in [g6[ the simple case of states with- 
out negative momenta so that there is no need to consider 
the energy degeneracy in the following expressions. The 
energy Fourier transform of ip(x — 0, t) is given by 

/+oo 
i>{ x = o, ty Et/h dt 
-00 

= h 1/2 (x = 0\E)(E\^(t = 0)}. 
In terms of r\ we can write 



(t) = -foj dE v{E)* nE - 

where fj(E) = r}{E)/ {j dt\ip{0,t)\ 2 f/ 2 ). Thus a "time 
operator" may be again identified with —ihdE, but now 
in the space of functions fj(E). This operator is different 
from Tab- 

In terms of the position and momentum operators the 
average presence time becomes 

,-v_ m ^(t^0)\p' 2 x + xp- 2 \yj{t = 0)) 

[1 2 ^( t = 0)\p-^(t = 0)) ' U 

provided the integrals exist. This is to be contrasted with 
the average time that can be defined using the current 
density J(x = 0, t), 

f°° dtJ(x = 0,t)t 

(*>j = r~ z r, r-rrr - m = OWabM* = o)) • 



F oo dtJ(x = 0,t) 



(2) 
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The different formal results are to be expected since, in 
the classical limit, ([!]) is associated with an average pres- 
ence time and (|^) with an average passage time. Note 
however that in both cases we can relate these mea- 
surements with the formal operator —iUBe, as has been 
pointed out all along. 

As stated above, no matter how one constructs general 
time operators, one is inevitably led to —iUds (plus, pos- 
sibly, a function of E) , but the physical interpretation of 
those operators, self-adjoint or associated with a POVM, 
is quite another issue. 



V. DISCUSSION 

In the sections above we have made reference to results 
and proposals of other workers in the subject. In this sec- 
tion, however, we intend to restate some of those results 
and discussions in the light of our description. To start 
with, none better than the initial paper of Aharonov and 
Bohm 0, where they introduce the operator Tab- They 
were not particularly concerned with the self-adjointncss 
or otherwise of the operator, but they did signal the prob- 
lem of the singularity for p — 0, which they brushed aside 
by stating that all components of the interesting wave- 
functions would be of high enough momentum. Insofar 
as that is the case, then there is indeed no problem in 
obtaining a probability distribution for measurements of 
the operator. Their real point of interest in introducing 
Tab was to show that there is no reason inherent in the 
principles of quantum mechanics for the energy of an ob- 
served system not to be observed in as short a time as 
one pleases, where the time is measured in an observer 
system. The topic of measurement limitations is retaken 
in much later work of Aharonov and others [ ^7f , where 
they make the distinction between direct and indirect 
measurements of time of arrival, indirect measurements 
being associated with Tab (or with regularized variants), 
and direct measurements with couplings of the particle 
with other degrees of freedom that in the classical limit 
would measure the time of arrival. They describe a lower 
bound for the time uncertainty of direct measurements 
inversely proportional to the typical energy of the particle 
(not to the energy uncertainty). In fact the same depen- 
dence has been recently described, and traced back to 
the lower bound of the energy spectrum, for the indirect 
measurement case pgfl . But the time uncertainty refers 
to an ensemble and does not preclude a precise determi- 
nation of the time of arrival in an individual measure- 
ment, in other words, this does not involve a limitation 
of the quantum theory to handle the arrival time concept 
and prescribe intrinsic arrival time distributions (formu- 
lated without explicit recourse to a measuring apparatus 
or an additional "clock" degree of freedom). Numeri- 
cal results have shown essential agreement between the 
(intrinsic) POVM distribution and operational distribu- 
tions given by phcnomenological screen models, except 



in rather pathological cases p9| , p0[ . 

Going back in time again, we come to the work of 
Allcock As we have signalled above, he advocated 
the use of the energy representation. Even more, he sig- 
nalled the need for negative energies in a particular way. 
If we indeed want to measure arrival times, we would 
have to prepare wavepackets that are localized on one 
side of the point of arrival (x — xq, say) for all times 
before a given one, t < to. This directly entails that the 
temporal Fourier transform of the wavefunction will ex- 
hibit all positive and negative energies, save for a set of 
null measure. In a way, it could be argued that the con- 
ceptual setup for measuring time of arrival requires by 
itself the introduction of negative energies. Furthermore, 
the proof he proposes for the need for all positive and 
negative energies is also illuminating: it goes the way of 
the Paley-Wiener-Titchmarsh theorem and the unique- 
ness theorem for analytic continuations. This is precisely 
the same route taken to show the non-orthogonality of 
any set of functions that could be considered as (gener- 
alized) eigenfunctions of a time operator. In that proof 
it is clear that the introduction of negative energies also 
eliminates the problem of non-orthogonality. However, 
his interpretation of the non-orthogonality of (general- 
ized) eigenfunctions for the time of arrival as precluding 
any measurement thereof is at odds with the operational 
approach within which the POVM idea takes place. 

A number of more mathematical papers has also been 
published, with a more functional analytical approach 
than those mentioned before |p|, ^0| , |l4| . 

The axiomatic approach of Kijowski Jll| has been men- 
tioned all along, since its guiding idea, i.e., characterizing 
the physical content of a concept such as time of arrival in 
a probability distribution, is crucial for our purposes, and 
clears the path for the introduction of the POVM as the 
physically relevant object. Kijowski also emphasized the 
usefulness of the energy representation, and decomposed 
the Hilbert space of the free particle as the direct sum of 
the subspaces of positive and negative momenta. He con- 
structs both the operators T± (in his notation, f^, if Q 
is the point X3 — 0) by comparison with the probability 
distribution, as we do. He then points out that T + T_ 
admits no self-adjoint extension, and goes no further that 
way. However, he also introduces a particularly interest- 
ing object, retaken more recently by Delgado and Muga 

l| (see also @), namely f' = f + ® (-^-)- He then 
asserts that this operator is essentially self-adjoint. Let 
us examine this question from the point of view taken in 
this paper. As we have signalled, both T + and T_ are iso- 
morphic to the momentum operator on the half-line, and 
their deficiency indices are (1, 0) for both of them. There- 
fore, — T- is also a maximally symmetric operator, with 
deficiency indices (0, 1). The total operator T + T^j , 

symmetric on a dense subspace, has deficiency indices 
(1,1), so it is not self-adjoint, although it admits self- 
adjoint extensions. Contrary to what might be inferred 
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from some claims in 12 1, there are infinite self-adjoint 
extensions, parametrized by a G [0,27r). The extension 
T' a is defined on D(T' a ), which is the set of elements of 
H+ 8 H- of the form + Xe- E / h , tp^ + \e la e~ E ' h ), 

where ip± £ D(T±), and A £ C. The action on its do- 
main is simply that of the couple of differential operators 
(—ihdE,i?idE)- The boundary conditions, however, are 
different. In terms of the momentum representation, the 
functions in the domain of T' a can have p 1 ^ 2 behaviour, 
for instance. At any rate, the probability distribution de- 
duced from these operators is not covariant under time 
shifts, so their physical interpretation is rather clouded 
unless the state belongs to only one of the two subspaces 
(positive or negative momenta) p0| . 

Another interesting attempt to construct a self-adjoint 
operator related to the time of arrival is that of Grot, 
Rovelli and Tate p6[ . In essence, they give a realiza- 
tion of the intuition of Aharonov and Bohm jjj that for 
high momenta the p — singularity of the differential 
operator is irrelevant, a point later taken up by Paul Q 
as well. However, this generalized self-adjoint operator 
has also some drawbacks from the point of view of its 
interpretation pO||r^ ]. From the point of view taken in 
this paper, the problem with the proposal of Grot, Rov- 
elli and Tate is that the strong operator limit of their 
deformed operator does not exist when their regulator is 
made to disappear (this can be seen, for example, com- 
puting the norm of the action by the regulated operator 
on any gaussian state, and then trying to eliminate the 
regulator - the norm blows up). Therefore, the results ob- 
tained with the regulator cannot be at all extrapolated 
to the unregulated case. 

Closest to the work carried out in this paper is the 
presentation by Busch, Grabowski and Lahti |||(| . How- 
ever, they refer the idea of using Naimark's dilation in 
this context to the paper of Werner |b| , where he con- 
structs "screen observables" using a version of imprimi- 
tivity adequate for the case of POVMs. Our analysis is 
much more simplified and explicit. 

In conclusion, we have given a description of the 
POVMs associated with the momentum operator on the 
half-line and the time operator of Aharonov and Bohm in 
terms of Naimark's dilations. More concretely, in terms 
of what for the momentum operator is the natural dila- 
tion. We have clarified several issues that have been the 
matter of debate or confusion in the literature as of late, 
in particular the deficiency indices and constructions of 
self-adjoint extensions of Tab, and we have rederived the 
probability density for times-of-arrival derived axiomati- 
cally by Kijowski. We have also made a proposal for the 
construction of the probability distributions of a wide 
class of time operators, for which it has been shown that 
they are not in general related to the arrival time, in spite 
of the formal agreement with the arrival time operator of 
the free motion case. 
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